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A UNIFIED  APPROACH  FOR  MODELING 
INELASTIC  BEHAVIOR  OF  STRUCTURAL 
METALS  UNDER  COMPLEX  CYCLIC 
LOADINGS 


1 INTRODUCTION 


Background 

The  importance  of  analyzing  seismic  structures  for 
inelastic  response  and  low-cycle  fatigue  lias  been  identi- 
fied by  several  investigators.1  b These  structures  are  de- 
signed to  respond  inelastieally  during  strong-motion 
earthquakes  both  for  economic  reasons,  and  to  provide 
effective  energy  dissipation  and  damping  through  plas- 
tic hysteresis.  Analysis  of  the  inelastic  response  of  such 
structures  is  complicated  by  the  cyclic  nature  of  earth- 
quake loads.  Current  structural  analysis  procedures  use 
material  constitutive  equations  which  are  only  adequate 
for  monotonic  loadings.  Under  cyclic  loading  condi- 
tions, inelastic  material  response  is  complicated  by  two 
factors:  ( 1 ) noncoincidence  of  tensile  and  compressive 
stress-strain  paths  resulting  in  hysteresis,  and  (2)  strong 
dependence  of  stress-strain  paths  on  prior  deformation 
history.  Tins  concept  is  illustrated  in  Figure  I,  which 
represents  a typical  stress-strain  response  of  a material 
subjected  to  a cyclic  deformation  history.  There  is  no 
unique  relationship  between  stress  and  strain  or  tan- 
gent modulus  and  stress  (or  strain);  these  quantities  are 
functions  of  their  time  histories.  Thus,  analysis  of  cy- 
clic inelastic  response  requires  a material  model  which 
can  adequately  simulate  the  observed  material  behavior 
and  its  prior  history  dependence. 

Hi.  V.  Berg  and  J.  L.  Strathu.  Anchorage  and  Alaska  Earth- 
quake (American  Iron  and  Steel  Institute,  1964). 

2S.  A.  Mahin  and  V.  V.  Bertero,  “Nonlinear  Seismic  Re- 
sponse 1 valuation  Charaima  Building."  Proceedings  ASCE. 
Journal  of  Structural  Division  (June  1974),  pp  1225-1242. 

3I  P Popov.  “Low  Cycle  Fatigue  of  Connections  and  De- 
tails." Proceed  inns  ASCE-IABSE  Joint  International  Confer 
cnee  on  Planning  and  Design  of  Tall  Buildings.  State-of-the-Art 
Report  No.  J.  Tech.  Comm.  No.  IX.  Vol  II  ( 1 972), pp  74 1-755. 

4(  \V.  Pinkliam.  “Procedures  and  Criteria  for  Earthquake 

Resistant  Design  Part  \."  Building  Practices  for  Disaster  Mitiga- 
tion. Building  Science  Scries  46  (National  Bureau  of  Standards. 
1973),  pp  I8S-208. 

■Mi  Bresler.  “Behavior  of  Structural  Elements  A Review,” 
Building  Practices  Jdr  Disaster  Mitigation.  Building  Science  Ser- 
ies 46  (National  Bureau  of  Standards,  1 973),  pp  286-35 1 . 

6N.  M.  Newmark  and  f . Rosenblueth,  Fundamentals  of 
Fartlu/uake  Fngineering  (Prentice  Mali.  Inc..  1971). 


Objective 

The  objective  of  this  report  is  to  develop  a unified 
approach  for  modeling  the  uniaxial  inelastic  behavior 
of  structural  metals  appropriate  for  irregular  cyclic 
loading  conditions.  Development  of  appropriate  materi- 
al parameters  required  for  such  a model  is  also  included 
in  this  objective.  The  model  will  be  digital-computer- 
oriented  and  suitable  for  numerical  procedures  in  struc- 
tural mechanics  and  failure  analysis. 

Approach 

It  is  convenient  to  distinguish  two  features  of  cyclic- 
inelastic  behavior  of  structural  metals  for  the  purpose 
of  modeling.  These  features  are  termed  the  “memory" 
phenomenon  (Chapter  2)  and  "plastic  hysteresis"  phe- 
nomena (Chapter  3).  Excellent  models  of  simulating 
memory  are  available7-!  1 and  the  one  best  suited  for 
the  present  purpose  will  be  adopted.  By  extending  re- 
cent studies'  2.  1 3 on  the  plastic  hysteresis  phenomena, 
a unified  approach  for  their  quantitative  description 
will  be  developed  and  combined  with  the  memory  mo- 
del to  obtain  a complete  characterization  of  material 
behavior.  Chapter  4 discusses  the  taterial  properties 
required  by  the  model  and  the  te;  is  for  determining 
those  properties.  Chapter  5 discus  es  the  synthesis  of 
memory  model  and  plastic  hystere  as  phenomena,  some 
appropriate  simplifications  in  the  model,  and  typical 
simulations  for  two  materials.  Chapter  6 outlines  three 


7W.  D.  Iwan.  "On  a Class  of  Models  for  the  Yielding  Beha- 
vior of  Continuous  and  Composite  Systems.”  Transactions  o] 
ASME,  Journal  of  Applied  Mechanics  (September  1967). 

8P.  C.  Jennings.  "Earthquake  Response  of  a Yielding  Struc- 
ture,” Proceedings  ASCE,  Journal  of  Engineering  Mechanics 
Division.  Paper  No.  4435,  EM  4 (August  1965),  pp  41-68. 

9J.  F.  Martin,  T.  H.  Topper,  and  11.  M.  Sinclair.  “Computer 
Based  Simulation  of  Cyclic  Stress-Strain  Behavior  With  Appli- 
cations to  f atigue, " Materials,  Research  and  Standards.  MTRSA  . 
Vol  1 1.  No.  2,  (February  1971),  pp  23-29. 

'°ll.  R.  Jhansalc  and  T.  11.  Topper.  “An  Engineering  Anal- 
ysis of  the  Inelastic  Stress  Response  of  a Structural  Metal  Un- 
der Variable  Cyclic  Strains,”  Cyclic  Stress-Strain  Behavior  t u 
alvsis.  Experimentation  and  Failure  Prediction . ASTM  S EP  5 1 9 
(American  Society  for  Testing  and  Materials  |ASTM) , 1973). 

llR.  M.  Wetzel,  A Method  of  Fatigue  Damage  Analysis. 
Technical  Report  No.  SR  71-107  (Ford  Motor  Company.  Aug- 
ust 1971). 

*2I1.  R.  Jhansalc,  “A  New  Parameter  for  the  Hysleretic 
Stress-Strain  Behavior  of  Metals,”  Transactions  of  ASME.  Jour 
nal  of  Engineering  Materials  and  Technology . Vol  97,  Series  It. 
No.  1 (January  1975),  pp  33-38. 

,3H.  R.  Jhansalc.  “A  Friction  Stress  Method  tor  the  Cyclic 
Inelastic  Behavior  of  Metals,”  Transactions  oj  the  Jrd  Interna 
tional  Conference  on  Structural  Mechanics  in  Reactor  Tech 
nology,  Vol  5,  t.5/4  (September  1975). 
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applications  of  structural  mechanics  and  failure  analy- 
sis for  this  research. 


STRAIN -TIME  INPUT 
CT 


STRESS- STRAIN  RESPONSE 


Figure  1.  Typical  cyclic  inelastic  material  response. 


Mode  of  Technology  Transfer 

The  results  of  this  study  will  impact  TM  5-80‘MO, 
Seismic  Design  for  Buildings . 1 •* 

2 MEMORY  PHENOMENON 

The  memory  phenomenon  can  be  best  described  by 
two  simple  examples,  as  shown  in  Figure  2.  If  the  ma- 
terial is  loaded  through  OA,  as  shown  in  Figure  2a.  then 
unloaded  through  AB.  and  finally  loaded  in  the  tensile 
direction,  it  apparently  "remembers"  the  previous 
point.  A.  of  unloading  and  follows  the  path  AC.  which 
would  be  the  path  if  the  material  were  not  unloaded  at 
A.  Thus,  there  is  a discontinuity  in  the  stress-strain  re- 
sponse during  the  loading  BC  at  A.  Similarly,  in  the 
second  example  (Figure  2b).  the  final  stress-strain  path, 
DF,  is  composed  of  two  discontinuities.  C and  A.  Seg- 
ments CA  and  AF  are  continuations  of  BC  and  OA.  re- 
spectively. The  memory  phenomenon  accounts  for  the 
internal  stress  distribution  caused  by  the  previous  load- 
ing history. 

Model  for  Memory 

A rheological  model  consisting  of  linear  elastic 
springs  and  solid  friction  sliders  such  as  the  one  shown 
in  Figure  3 is  commonly  used  to  simulate  the  memory 
phenomenon.  By  choosing  an  appropriate  number  of 
spring-slider  elements,  stiffnesses  for  the  springs,  and 
friction  stresses  for  the  sliders,  the  nonlinear  stress- 
strain  characteristic  of  the  material  can  be  matched  as 
accurately  as  desired  with  a piecewise  linear  approxi- 
mation. The  model  also  exhibits  the  so-called  Baushin- 
ger  effect,  which  is  characterized  by  a decrease  in  the 
yield  stress  level  in  the  reverse  direction  after  plastic 
deformation.  Based  on  this  rheological  model,  several 
formulations  for  simulating  the  memory  phenomenon 
have  been  proposed.1  5 • 1 <’  The  approach  which  uses  a 

1 ^Seismic  Design  for  Buildings,  TM5-809-10  (Department 
of  the  Army.  197.1).  This  I'M  is  also  published  by  the  Navy  and 
Air  Force  as  NAVI' AC'  P-355  and  Chapter  1 3 of  AFM88-3,  re- 
spectively. 

15J.  F.  Martin.  T.  It.  Topper,  and  CT  M.  Sinclair.  “Compu- 
ter Based  Simulation  of  Cyclic  Stress-Strain  Behavior  with  Ap- 
plications to  Fatigue,”  Materials,  Research  and  Standards. 
MTRSA , Vol  1 1 . No.  2,  (February  1971 ).  pp  23-2 9. 

16|<  M.  Wetzel.  A Method  of  Fatigue  Damage  Analysis, 
Technical  Report  No  SR  71-107  (l  ord  Motor  Company,  Aug- 
ust 1971). 
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(a)  <b) 

Figure  2.  Two  examples  of  memory  phenomenon. 


CT1  <Jn 


E 

~ 'AA/V  Elastic  Spring  with  Modulus, E 

■ w Solid  Friction  Slider,  with 
Friction  Stress,  <f 

ffi-i  < aj  < aj+1 

Figure  3.  Rheological  model. 

sel  of  simple  rules  is  ideally  suited  for  digital  computa- 
tion and  will  be  adopted.1  7 These  rules  will  be  illustrat- 
ed with  the  aid  of  the  following  example. 

Methods  of  Memory  Rules 

lhe  initial  stress-strain  curves  for  tension  and  com- 
pression ate  assumed  to  be  identical,  flic  stress-strain 


curves  can  be  approximated  by  a set  of  piecewise  linear 
segments  whose  lengths,  slopes,  and  number  are  suit- 
ably chosen.  Figure  4 illustrates  a live-segment  curve. 
The  memory  rules  are  as  follows: 

1 . The  availability  of  a segment  during  the  current 
loading  path  which  is  either  tensile  or  compressive  de- 
pends on  the  prior  loading  history  and  is  denoted  by  an 
“availability  coefficient"  in  that  direction. 

2.  The  absolute  sum  of  the  availability  coefficients 
in  the  tensile  and  the  compressive  directions  for  each 
segment  is  always  equal  to  two.  The  initial  availability 
coefficient  in  tension  or  in  compression  for  each  seg- 
ment is  one. 

3.  During  a specific  loading,  the  stress-strain  path  is 
defined  by  segments,  starting  with  the  first,  and  pro- 
ceeding in  consecutive  order  to  the  extent  that  the  seg- 
ments are  available  in  the  direction  of  loading  until  the 
desired  stress  or  strain  limit  is  reached. 

4.  The  availability  of  a segment  in  a given  direction 
(tension  or  compression)  decreases  to  the  extent  it  is 
used  in  that  direction,  but  increases  by  that  amount  in 
the  opposite  direction  in  conformity  with  tide  2. 

In  brief,  the  procedure  amounts  to  a simple  book 
keeping  operation  of  the  availability  coefficients  as  il- 
hist  rated  in  Figure  4.  When  the  tensile  loading  sequence 


o 
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Availability  Coefficient  of  Elements 


Loading 

Path 


c 

T 

c 

1.0 

1.  0 

1.  0 

2.0 

0.  2 

1.  8 

2.  0 

0.  2 

1.  8 

2.  0 

0.  2 

1.  8 

© 

T 

c 

1.0 

1.  0 

1.  0 

1.  0 

2. 0 0. 0 1. 3 0.  7 0. 0 2.  0 0.  2 


. 0 0. 0 2. 0 0.  0 2. 0 0.  0 2. 0 0. 0 2. 0 


1.  0 

I.  0 

0.0 

2.  0 

(S  - Start,  1'.  - End,  T - Tension.  C - Compression,  Unused  Elements  Underlined) 


Example  for  illustrating  memory  rules.  Reprinted  with  permission  of  Ford  Motor  Company  front  R.  M 
Wetzel.  .1  Mctfuhl  of  Fatigue  Damage  Analysis,  Technical  Report  No.  SR7I-107  (August  I ‘>71 ). 


ACT 


Ac-  stress  range 

AC  - STRAIN  RANGE 
Ca  -MEAN  STRESS 


(a)  Stable  closed  loop  behavior  under  mean  stress.  (b)  Identically  shaped  hysteresis  branches  (yield  range 

is  constant). 

Figure  5.  Masing  and  stable  responses  of  memory  model. 


OA  begins,  one  unit  of  each  segment  is  available  in 
tension  and  compression.  Hence,  the  stress-strain  path 
OA  consists  of  one-unit  lengths  of  elements  1 to  3 and 
0.8  unit  of  4;  at  this  stage,  point  A is  reached.  The  avail- 
ability coefficients  of  these  segments  in  tension  are 
reduced  by  the  amounts  of  these  element  lengths,  but 
at  the  same  time,  their  availability  coefficients  in 
compression  are  correspondingly  increased  by  the  same 
amount.  The  availability  coefficients  of  segment  5 are 
unaltered,  since  that  segment  was  not  used.  The  next 
loading  sequence.  AB.  which  is  compressive,  uses  the 
two  available  units  of  segment  1 and  1.3  units  of  seg- 
ment 2;  at  this  stage,  the  desired  limit  B is  reached.  The 
next  sequence,  BC.  which  is  tensile,  uses  the  two  units 
of  1 , 1.3  units  of  2,  none  of  3, 0.2  units  of  4,  and  one 
unit  of  5;  at  this  stage,  point  ('  is  reached.  Thus,  the 
memory  phenomenon  has  been  successfully  simulated 
by  recognizing  the  unloading  point  A and  producing 
stress-strain  path  AC.  which  is  the  continuation  of  OA. 

Memory  Model  Vs.  Real  Materials 

Under  constant  stress  or  strain  range  cycling,  the  rhe- 
ological model  (or  the  method  of  rules)  produces  fully 
closed  Its  steresis  stress-strain  loops,  even  in  the  presence 
of  a mean  stiess  or  mean  strain  as  illustrated  in  Figure 
5.  Although  real  materials  eventually  approach  this  sta- 


ble (or  saturation)  state  in  the  limit,  they  usually  ex- 
hibit a prior  transient  behavior  (to  be  fully  described  in 
Chapter  3)  during  which  hysteresis  paths  gradually 
change  from  cycle  to  cycle.  This  transient  behavior  is 
particularly  pronounced  in  the  presence  of  a mean 
stress. 

The  hysteresis  curves  produced  by  the  rheological 
model  (or  the  method  of  rules)  are  geometrically  simi- 
lar to  the  monotonic  stress-strain  curve  but  magnified 
by  a scale  factor  of  two.  This  occurs  because  in  order 
to  reverse  the  sliding  direction  lot  a slider,  it  is  neces- 
sary to  applj  a stress  range  which  is  twice  that  ol  the 
slider’s  friction  stress.  This  property,  originally  postu- 
lated by  Masing.  is  often  referred  to  as  Masing  postula- 
tion (or  Masing  behavior). IS  It  also  implies  that  the  lin- 
ear elastic  portion  or  yield  range  (which  is  twice  the 
yield  strength)  is  independent  of  the  hysteresis  loop  size 
(or  strain  or  stress  amplitude).  For  example,  if  various 
sized  hysteresis  loops  produced  by  the  model  ate  super- 
imposed on  their  lower  tips,  the  upper  hysteresis  paths 
are  coincident,  as  shown  in  Figure  5b.  Many  real  mater- 


18(i  Masing,  "Figenspannungen  und  Veriest igung  bein 
Messing,”  Proceedings  of  the  2nd  International  < onyress  of  1 /> 
plied  Mechanics,  Zorich  IV26  f 1926),  pp  332-335. 


lals,  even  in  the  stabilized  state,  exhibit  a variation  of 
yield  range  as  a function  of  loop  st/e.  These  materials 
are  denoted  as  non-Masing  type. 

These  two  features  (the  transient  and  the  non-Masing 
characteristics  which  distinguish  real  materials  from 
the  memory  model)  are  designated  as  plastic  hysteresis 
phenomena  (see  Chapter  3). 


3 PLASTIC  HYSTERESIS  PHENOMENA 

A recent  study  of  several  structural  metals19  showed 
that  the  various  cycle-dependent  transient  phenomena 
and  non-Masing  behavior  can  be  characterized  in  terms 
of  a single  history-dependent  stress  parameter.  The 
present  study  has  yielded  additional  insight  into  the 
quantitative  aspects  of  this  parameter  and  its  history 
dependence.  This  chapter  describes  these  findings  and 
provides  a suitable  model  for  characterizing  the  plastic 
hysteresis  phenomena. 

Observed  Transient  Phenomena 

The  transient  phenomena  are  classified  into  four 
types  (Figure  6)  and  are  designated  as  cyclic  hardening, 
cyclic  softening,  cyclic  relaxation,  and  cyclic  creep.  As 
illustrated  in  Figures  6a  and  6b,  cyclic  hardening  or 
softening  is  denoted  by  a cycle-dependent  increase  or 
decrease  in  stress  range  under  constant  strain  range  cy- 
cling. Similarly,  under  constant  stress  range  cycling, 
these  two  phenomena  manifest  as  changes  in  the  strain 
range.  Cyclic  relaxation  is  denoted  by  a cyclic-depen- 
dent decrease  in  the  absolute  mean  stress  under  con- 
stant strain  range  cycling,  whereas  cyclic  creep  is  denot- 
ed by  a cycle-dependent  increase  in  the  absolute  mean 
strain  under  constant  stress  range  cycling.  Figures  6c 
and  6d  schematically  illustrate  cyclic  relaxation  of  a 
tensile  mean  stress  and  cyclic  creep  in  tension,  respec- 
tively. A material  can  simultaneously  exhibit  cyclic 
hardening  or  softening  and  cyclic  creep  or  relaxation. 

Extensive  experimental  observations  of  the  transient 
behavior  of  many  different  structural  metals  made  by 


19|l.  R.  J ha  male.  “A  New  Parameter  fur  the  Hysteretic 
Stress-Strain  Behavior  of  Metals,"  Transactions  of  ASM  l-\  Jour- 
nal oj  Engineering  Materials  ami  Technology , Vol  97.  Series  li, 
No  I (January  1975).  pp  33-38. 


several  investigators  during  the  past  two  decades-11--1 
can  be  generally  summarized  as  follows.  Soft  metals 
(e.g..  hotworked)  cyclically  harden  and  hard  metals 
(e.g.,  coldworked)  cyclically  soften.  However,  initially 
soft  but  strain-aged  materials,  such  as  hotworked  com- 
mon structural  steels,  exhibit  cyclic  softening  at  lower 
strain  amplitudes  and  cyclic  hardening  at  higher  strarn 
amplitudes.  This  apparent  anomalous  softening  is 
believed  to  be  caused  by  the  gradual  unpinning  of  dis- 
locations which  were  locked  by  interstitial  impurity 
atoms  during  the  strain-aging  process.  One  feature  of  a 
softening  phenomenon  that  should  be  of  concern  in 
deformation  analysis  is  the  fact  that  what  appears  to  be 
a linear  elastic  response  during  the  initial  period  of  cyclic 
loading  will  eventually  develop  into  inelastic  hysteresis 
response  with  continued  cycling  as  illustrated  in  Figure 
7.  Thus,  analysis  based  on  the  initial  stress-strain  char- 
acteristic could  result  in  nonconservative  estimates  of 
deformation  response. 

The  initial  cyclic  hardening  or  softening  phenomenon 
is  transient,  and  a state  of  saturation  is  approached  at  a 
continuously  decreasing  rate  with  cycles.  This  stable  or 
saturated  state  is  usually  maintained  with  continued  cy- 
cling at  constant  amplitude  until  fatigue  cracks  initiate. 
However,  each  time  the  stress  or  strain  amplitude  (or 
limit)  is  changed  from  one  level  to  another,  a transient 
behavior  ensues  even  though  the  material  has  been  in- 
itially saturated.  This  transient  behavior  is  comprised 
of  both  hardening  or  softening  and  relaxation  or  creep. 
Specific  data  illustrating  these  various  aspects  of  ob- 
served transient  behavior  are  discussed  in  the  following 
sections  dealing  with  the  quantitative  characterization 
of  these  phenomena. 

Analysis  of  Hysteresis  Loops 

A hysteresis  loop  has  two  branches  the  upper  branch 
corresponds  to  the  tensile  loading  direction  and  the 
lower  branch  corresponds  to  tin  ^impressive  dueclion. 
During  the  transient  conditions,  a hysteresis  loop  is  not 
closed,  and  the  hysteresis  biaim  s me  nm  identical 
(see  Figure  6);  however,  iii  the  satui  tied  state,  the  hys- 

I Coffin,  lr  ami  l i favernefe  *Th* tydh  Strain 
ing  and  Fatigue  ,,i  Metals  ' Tranaaetiom  oj  r*n  \l<  tallurgu  al 
Society.  I /.)//■  V.tl  71  5 t(ktohei  19S  o pp  ‘ M s -i 

R I uler  and  J Morrow  . Cede  Ihrfn-nJent  Strew  Strain 
Behavior  oj  Metals  I A \M  Report  N ‘ ' it  niversily  ol  till 
liois.  .Man'll  19h.t| 

22R  \V  Landgr.il . I i du  in  formation  ami  I attain  tijllarj 
cnctl  Steels  KV  \M  Report  J2(l  <1  niversily  ol  Illinois  l'*l,Xi 

23s  Kesltaven.  .Sour,  .Srin/l,  \ on  the  information  ami  hrai 
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Figure  7.  Gradual  development  of  inelastic  hysteresis  from  initially  elastic  behavior  in  cyclically  softening  materials. 


teresis  loop  is  fully  closed,  and  the  two  hysteresis 
branches  are  identical,  one  being  the  mirror  image  of 
the  other.  Each  branch  is  designated  as  a "reversal,”  and 
cyclic  progress  is  denoted  in  terms  of  the  number  of 
reversals.  The  first  reversal  is  the  hysteresis  branch  im- 
mediately following  the  first  monotonic  loading  at  the 
beginning  of  the  test.  As  a rule,  the  first  monotonic 
loading  path  is  treated  separately  from  the  other  cyclic 
loading  paths,  and  no  reversal  number  is  assigned  to  it. 

A recent  study  of  several  structural  metals-4  showed 
that  transient  changes  in  the  hysteresis  branches  are  es- 
sentially caused  by  changes  in  the  lengths  of  the  initial 
“linear"  elastic  parts.  The  slope  of  the  "linear”  elastic 
portions  and  the  shape  of  the  nonlinear  portions  remain 
virtually  unchanged  (see  Figure  8,  in  which  the  first 
three  hysteresis  branches  of  a typical  cyclic  hardening 
situation  are  compared).  The  change  in  the  elastic  part 
is  designated  as  "yield  range  increment"  (VRI),  since  it 
denotes  a change  in  the  yield  range  of  the  material.  The 
yield  range  is  twice  the  yield  strength  because  of  the 
scale  factor  of  two  between  the  monotonic  and  hyster- 
etic  stress-strain  paths  (see  Chapter  2).  Physically,  this 
change  in  the  yield  strength  (or  range)  is  caused  by  sub- 
structural  changes  associated  with  reversed  plastic 
straining  in  a cyclic  situation.  Thus,  transient  phenom- 
ena can  be  described  in  terms  of  a simple  parameter 
VRI  or  "yield  strength  increment"  ( V SI),  which  is  one- 
hall  of  YRI  I he  elastic  modulus  can  he  assumed  to  re- 


-4lf.  R ihansitle.  A New  Parameter  for  the  Hysteretic 
Stress  Strain  behavior  of  Metals.  Transactions  oj  ASME , Jour- 
nal ol  Ungineerinx  Materials  ami  Technology,  Vol  97,  Series  II. 
No  1 (January  1 975),  pp  33-38. 


main  unchanged  during  cyclic  loading.  The  invariant 
nonlinear  portion  of  the  hysteresis  branch  that  is  un- 
affected by  the  cyclic  plastic  straining  represents  an 
intrinsic  material  property.  Denoting  such  an  intrinsic 
stress-strain  curve  as  "skeleton  stress-strain  curve,”  the 
invariant  nonlinear  hysteresis  branch  will  be  geomet- 
rically similar  to  the  skeleton  stress-strain  curve,  but 
magnified  by  a scale  factor  of  two.  Chapter  4 discusses 
determination  of  the  skeleton  stress-strain  curve. 

A reversal  to  reversal  increase  or  decrease  in  YRI  de- 
notes hardening  or  softening,  respectively,  whereas  a 
relative  difference  in  YRI  between  the  upper  and  lower 
hysteresis  branches  describes  cyclic  creep  or  relaxation, 
depending  on  the  controlled  limits  (stress  or  strain),  as 
schematically  illustrated  in  Figure  9.  The  relative  dif- 
ference in  YRIs  between  two  consecutive  reversals  de- 
creases as  saturation  is  approached.  Therefore,  in  a sat- 
urated hysteresis  loop,  the  YRIs  in  the  upper  and  lower 
branches  are  equal  in  magnitude.  In  most  materials,  the 
saturation  value  of  YRI  is  dependent  on  the  stress  or 
strain  amplitude.  These  materials  were  designated  as 
non-Masing  type  in  Chapter  2. 

As  part  of  the  present  investigation,  all  the  observa- 
tions discussed  above  have  been  verified  for  a large 
number  of  structural  metals  ranging  widely  in  strengths 
and  cyclic  characteristics.  Analysis  of  available  data  on 
the  elevated  temperature  combined  cyclic  and  hold 
period  (static  relaxation/creep)  tests  on  two  high-tem- 
perature alloy  steels  also  indicate  that  both  cyclic  tran- 
sient and  static  creep/relaxation  behavior  can  be  de- 
scribed in  terms  of  changes  in  YRI.  Table  1 summarizes 
the  materials  studied  and  observations  about  them.  Only 
typical  results  are  discussed  in  the  following  section. 
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(a)  Three  hysteresis  reversals  of  (b)  Three  reversals  are  superim-  (c)Reversals  1 and  2 are  translated 

cyclic  hardening.  posed  on  their  initial  points  Aj  along  linear  elastic  slope  to 

Bj.  and  Cj,  respectively.  match  nonlinear  portions.  Curve 

AjXDCB  is  a doubled  skeleton 
stress-strain  curve;  linear  parts 
Aj.  Bj.  and  AjCj  are  YRIs  of 
reversals  2 and  3,  respectively. 
YR1  of  reversal  1 is /tero. 

Figure  8.  Comparison  of  first  three  hysteresis  branches  of  a typical  hardening  situation. 


(a)  Cyclic  hardening.  (b) Cyclic  relaxation.  (c)  Cyclic  creep. 

Figure  9.  Reversal  to  reversal  variation  in  YRI  during  transient  behavior.  AB.  CD,  and  IT  are  YRIs  corresponding 
to  reversals  AC,  CM,  and  MCI.  Portions  BC.  DM.  and  1(1  are  identical  in  shape  and  arc  described  by  dou- 
bled skeleton  stress-strain  curves. 


Table  1 

Summary  of  Hysteresis  Loop  Shape  Study  of  Structural  Metals 


C onditions  of  l est 


Initial  Yield  Exhibits  Cyclic  Types  of  Is  Material 


No. 

Material  Name 

Strength 
in  ksi  (MHa) 

Hardening  or 
Softening 

Temp. 

History 

Hysteresis 
Loops  Studied 

Masing  or  Non- 
Masing  Type? 

Data  Source 

STEELS 

1 

SAL  1018  Steel 
normalized 

33 

(228) 

ll/S 

RT 

Constant  T&St 

Strain 

Amplitude 

N 

a,  b 

2 

A STM  A-36 
steel 

40 

(276) 

H/S 

RT 

" 

r&si 

N 

c 

3 

SAL  1045  steel 
normalized  \ 

45 

(311) 

H/S 

RT 

” 

T&St 

N 

d 

4 

A STM  A -440 
steel  (Man-Ten) 

46 

(317) 

H/S 

RT 

** 

Si 

N 

e 

5 

lladfied  steel 

60 

(414) 

11 

RT 

St 

N 

f 

6 

LN25  steel 

Q&T  (450°C)  British 

90 

(621 ) 

S 

RT 

T 

N 

£ 

7 

SAL  1045  steel 

C2&T  ( 1 200°  1'  t649°(  |) 

92 

(635) 

s 

RT 

” 

T&T 

N 

li 

8 

SAL  4340  steel 

100 

(690) 

s 

RT 

" 

S( 

N 

i 

9 

Cast  8630  steel 

103 

(711) 

s 

RT 

liSt 

N 

J 

III 

RCX'-IOO  steel 

130 

(897) 

s 

RT 

St 

N 

k 

1 1 

SAL  1045  steel 
Q&T  720” L 1 382“C| 

185 

(1277) 

s 

RT 

T&St 

N 

h 

12 

Mar  200  steel 

215 

(1484) 

s 

RT 

St 

N 

h 

13 

SAL  1045  steel 

((3&  r (500"  1 1 260° C'l  1 

245 

(1691) 

s 

RT 

" 

T&St 

N 

h 

14 

SAL  1045  steel 
y&T(360°l  1 1 82°C|  > 

260 

(1794) 

s 

RT 

r&si 

M 

h 

15 

SAL  1045  steel 
as  quenched 

265 

(18291 

SI 

RT 

" 

St 

M 

d 

16 

Ausformed  III  1 steel 

295  r 

Si 

RT 

- 

St 

M 

1 

(2036) 

365C 

(2519) 


I 


! 
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Table  I (Continued) 

Summary  of  Hysteresis  Loop  Shape  Study  of  Structural  Metals 


Conditions  of  Test 


, 

No. 

Material  Name 

Initial  Yield 
Strength 
in  ksi  (MPa) 

Exhibits  Cyclic- 
Hardening  or 
Softening 

Types  of 
Hysteresis 

Temp.  History  Loops  Studied 

Is  Material 
Masing  or  Non- 
Masing  Type? 

Data  Source 

17 

A1SI  304 
Stainless  Steel 

10 

(69) 

II 

1 1001-'  constant  T&Sl 

strain  cy- 
cling & 
static  creep 
and  relaxation 
during  holds. 

N 

m 

IK 

2 1/4  CR-l  Mo  steel 

26 

(179) 

11 

9501-  " T&St 

N 

n 

OTHER  METALS 

Oi  l 1C  copper 
annealed 

03 

(16.7) 

H 

RT  constant  St 

strain  cycling 

N 

b 

20 

2024  T4  aluminum 

55 

(380) 

11 

RT  ” T&St 

M 

d.  o 

21 

j 

7075- 16  aluminum 

68 

(469) 

H 

RT  " T&St 

M 

d 

Legend:  II  Hardening;  S Softening;  I Transient;  St  Stable;  M Masing  type,  N Non-Masing  type;  RT  Room  Temperature. 


a S.  keshaven.  Some  Studies  on  the  Deformation  and  Fracture  of  Normalized  Mild  Steel  Under  Cyclic  Conditions.  PhD.  Disserta 
tion  (University  of  Waterloo.  December  1966). 

b H.  R.  Jhansaie.  Inelastic  Deformation  and  Fatigue  Response  of  Spectrum  Loaded  Strain  Controlled  Axial  and  Flexural  Members. 
Ph  D.  Dissertation  (University  of  Waterloo,  March  1971). 

c J I*  Martin.  Cyclic  Mechanical  Tests  and  an  Appropriate  Analytical  Stress-Strain  Model  for  A -36  Steel . Technical  Report  M-36/ 
AD780H02  (U.  S.  Army  Construction  engineering  Research  Laboratory,  May  1974). 
d H.  R.  Jhansaie.  unpublished  experimental  studies,  11.  I\  Moore  fracture  Lab  (University  of  Illinois.  1973-74). 

e P.  C.  Rosenberger.  Fatigue  Behavior  of  Smooth  and  Notched  Specimens  of  Man-  Ten  Steel.  M S Thesis  (T&  AM  Department. 

University  of  Illinois,  1968). 

f R C.  Rice  and  R.  1.  Stephens,  Overload  Effects  on  Sub-Critical  Crack  Growth  in  Austenitic  Manganese  Steel,  ASTM  STP  536 
(1973). 

g P.  W.  J.  Oldroyd,  D.  J.  burns,  and  P.  P.  Benhain.  “Strain  Hardening  and  Softening  of  Metals  Produced  by  Cycles  of  Plastic  De- 
formation,” Proc.  Institution  of  Mechanical  Engineers,  Vol  150  (1965-66). 
h R.  W.  Landgraf.  unpublished  experimental  data, Graduate  Research  Program  (T&AM  Department.  University  of  Illinois.  1966-68) 

i R W.  Smith,  M.  H.  Hirschberg,  and  S.  S.  Man  son.  Fatigue  Behavior  of  Materials  in  Low  and  Intermediate  Life  Range,  NASA 

Technical  Note  D-1574  (National  Aeronautics  and  Space  Administration.  April  1963). 
j D.  1 Ditt mer  and  M.  R Mitchell.  Material  Characterization  of  Cast  S6J0  Steel  Monotonic  and  Cyclic  Stress-Strain  Behavior 
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o JoDean  Morrow,  unpublished  data  from  graduate  course  notes  ( f&AM  Department.  University  of  Illinois.  Spring  1973) 
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Observed  Typical  Yield  Range 
Increment  (YRt)  Behavior 

Most  materials  can  be  classified  into  three  general 
groups,  depending  on  the  type  of  initial  transient  behav- 
ior they  exhibit.  These  ate  ( 1 ) materials  which  exhibit 
initial  cyclic  hardening,  (2)  materials  which  exhibit  ini- 
tial cyclic  softening,  and  (3)  materials  which  exhibit  in- 
itial cyclic  softening  or  hardening  depending  on  the 
strain  amplitude.  The  observed  YR1  behavior  will  be 
discussed  in  terms  of  these  three  groups. 

First  Reversal  and  Saturation  YRIs 

The  values  of  YRI  corresponding  to  the  first  reversal 
and  cyclic  saturation  state  are  plotted  as  functions  of 
strain  amplitude  for  three  typical  materials  (each  rep- 
resenting a different  group)  in  Figures  10,  11.  and  12. 
l ot  each  material,  the  YRI  corresponding  to  the  state 
which  has  the  minimum  linear  elastic  portion  in  the 
hysteresis  curve  is  assigned  a /.ero  value.  The  YRIs  of 
all  other  hysteresis  branches  are  measured  with  respect 
to  this  datum  and  therefore  are  always  positive. 

In  the  case  of  2024-14  aluminum,  which  is  initially 
soft  and  exhibits  cyclic  hardening,  the  saturation  YRI 
corresponding  to  the  hardened  slate  is  larger  than  the 
first  reversal  value  shown  in  Figure  It).  The  first  rever- 
sal value  is  approximately  /ero.  and  both  the  initial  and 
saturation  values  of  YRI  are  essentially  independent  of 
strain  amplitude.  Consequently  , the  hysteresis  branches 
of  various  sized  saturation  loops  are  identical  in  shape, 
as  shown  in  Figure  13,  in  which  fully  reversed  hystere- 
sis loops  are  superimposed  on  their  lower  tips.  As 
shown  in  Table  1 . 2024-T4  is  one  of  the  few  materials 
observed  that  exhibit  Masing  behavior.  A comparison 
of  the  saturation  YRIs  obtained  under  constant  ampli- 
tude cycling  and  mixed  print  cyclic  history  shown  in 
Figure  10  indicates  that  the  saturation  YRI  is  relatively 
insensitive  to  prior  cyclic  history. 

In  the  case  of  SAI  1045  quenched  and  tempered 
steel,  which  is  initially  hard  and  exhibits  cyclic  soften 
mg,  both  the  first  reversal  and  saturation  YRIs  are  de- 
pendent on  strain  amplitude,  as  illustrated  m I igure  1 1 
Consistent  with  this  softening  phenomenon,  the  sat 
illation  YRI  ,s  less  than  the  lust  reversal  N Rl  ( old- 
worked  structural  steels  are  generally  expected  to  ex 
Inbit  these  characteristics 

In  the  case  of  common  structural  steels  tiniliallv 
hot  worked)  which  exhibit  initial  cyclic  softening  or 
hardening  depending  on  strain  amplitude  the  first  ic 
versa!  and  saturation  YRIs  vary  with  slum,  as  tvpically 
shown  for  A-36  steel  in  Figure  12  The  initial  high  val 


tie  of  YRI  corresponding  to  small  strains  is  due  to  the 
strain-aged  condition  of  the  material.  The  sudden  drop 
in  the  initial  YRI  on  the  application  ol  a small  mono- 
tome  strain  is  due  to  the  unpinning  of  dislocations 
which  were  locked  by  interstitial  impurity  atoms  dur- 
ing strain  aging.  This  phenomenon  is  also  manifested 
in  the  "upper  and  lower  yield  point”  behavioi  ol 
steels.  At  smaller  strain  levels,  where  the  first  reversal 
YRI  is  larger  than  the  saturation  YRI,  the  material  ex- 
hibits cyclic  softening  and  higltei  stiam  levels  when* 
the  saturation  YRI  is  larger  than  the  first  reveisal.  T RIs 
are  associated  with  cyclic  hardening.  The  non  Masing 
behavior,  signified  by  the  dependence  of  salutation 
YRI  on  strain  amplitude,  is  well  illustrated  in  Figure 
14.  in  which  saturated,  fully  icversed  hyxteiexis  loops 
ot  various  sizes  are  compareu  by  superimposing  them 
on  their  lower  tips  The  upper  branches  do  not  match: 
however,  by  sliding  the  individual  loops  along  the  elas 
tic  slope,  the  upper  branches  can  be  matched  extreme- 
ly well,  as  shown  in  Figure  15.  The  difference  in  the 
linear  elastic  portions  between  different  sized  loops  is 
due  to  the  dependence  of  YRI  on  strain  amplitude. 
As  shown  in  Figure  12.  the  saturation  YRI  appears  to 
be  somewhat  sensitive  to  prior  cyclic  history  at  smaller 
strains  and  not  at  higher  strains.  Past  inicrostructural 
stud  res  25  have  shown  that  loss  and  high-strain  levels 
are  generally  associated  with  widely  different  substiuc- 
tures;  it  is  rather  difficult  to  revert  from  the  one  corres- 
ponding to  the  higher  strain  level  to  the  one  cot  res- 
ponding to  lower  levels,  when  higher  strain  levels  art 
preceded  by  lower  one,  but  not  vice  versa. 

Variation  of  YRI  During 
Cyclic  llardcning/Si  iftcning 

To  measure  the  variation  ol  YRI  during  cyclic  hat- 
dening  softening,  it  is  necessary  to  analyze  the  tiansi 
cut  hysteresis  loops  under  constant  amplitude  cycling 
Such  transient  hysteresis  loop  data  were  available  for 
three  materials  representing  the  three  general  groups 
mentioned  previously  SAL  1045  normalized  steel. 
2024  14  aluminum,  and  SAI  1045  quenched  and  tem- 
pered (720°F  |37S°t'|  ) steel.  It  was  observed  that  in 
all  three  materials,  the  rate  of  YRI  change  pet  reveisal 
at  a given  stiam  amplitude  decreased  as  saturation  was 
approached  and  the  number  ol  reveisals  increased. 
Therefore,  the  following  general  expression  lot  the  rate 
ot  change  of  yield  strength  was  assumed  for  analyzing 
the  data 

“5t  t teller  amt  < l,imi.  "(‘ychc  Stress  Strain  Response 
ol  Ki  t Metals  .nut  Alloys.  I Phenomenological  I xperimenlv. 
II  Dislocation  Structures  amt  Mechanisms."  I<  la  llctallursica. 
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Figure  10.  First  reversal  and  saturation  YRIs  of  2024-T4  aluminum  expressed  as  functions  of  strain  amplitude 
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I nsi  reversal  and  saturation  YRIs  ol  SAF  1 045  quenched  and  tempcied(720  I |378°C|  ) steel  expressed 
as  functions  of  strain  and  amplitude 
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where  = rate  of  change  of  yield  strength  per  re- 

versal 


oy  = current  yield  strength 


°ys  = saturation  value  of  yield  strength 
R = current  reversal  number 


Quantitative  Characterization 
of  Plastic  Hysteresis  Phenomena 

To  obtain  a complete  characterization  of  the  various 
transient  phenomena,  one  must  know  the  variation  of 
YR1  (or  yield  strength)  both  from  reversal  to  reversal 
and  within  each  reversal.  The  discussion  presented  in 
the  previous  section  provides  a basis  for  the  former  var- 
iation: however,  it  is  relatively  difficult  to  measure  the 
variation  of  YRI  within  a given  reversal,  and  no  such 
experimental  data  are  available.  Therefore,  a physically 
motivated  intuitive  model  will  be  assumed  for  the  lat- 
ter variation. 


C||s  = coefficient  of  hardening  or  softening 
and  in  and  p are  suitable  exponents. 

Eq  1 applies  to  a constant  amplitude  cycling  situa- 
tion. and  the  values  of  Chs.  m,  and  p are  generally 
functions  of  strain  amplitude.  If  the  equation  is  expect- 
ed to  yield  both  the  magnitude  and  sign  of  the  rate  of 
yield  strength  change,  then  possible  values  of  m will  be 
restricted  to  odd  integers.  If  only  magnitude  is  obtained, 
then  such  a restriction  of  the  possible  value  of  m is  not 
necessary.  By  integrating  the  above  equation,  an  equa- 
tion relating  yield  strengths  and  number  of  reversals 
can  be  obtained  as  in  Eq  2: 


To  examine  how  well  Eq  2 represents  the  observed 
behavior  for  a specified  set  of  m and  p,  it  is  only  neces- 
sary to  evaluate  the  left-  and  right-hand  expressions  at 
various  reversals  and  to  examine  the  linearity  of  their 
relationship  on  an  x-y  plot.  Four  values  of  m and  p (0, 
1/2.  1,  and  2)  were  considered,  which  resulted  in  16 
combinations  for  examination.  The  best  combination 
varied  from  material  to  material  and  between  strain  lev- 
els for  each  material;  however,  the  one  combination 
which  provided  a consistent  and  reasonably  good  repre- 
sentation between  the  three  materials  and  different 
strain  levels  was  m = 1/2  or  1 , and  p = 1 . Figures  16, 17, 
and  18  are  typical  plots  of  expressions  of  the  equation 
for  m = 1 , p = 1 for  the  three  materials. 

From  the  above  analysis,  it  is  reasonable  to  assume 
that  under  constant  strain  amplitude  cycling,  the  change 
rate  of  yield  strength  is  proportional  to  the  difference 
between  the  saturation  and  current  values  of  yield 
strength,  and  inversely  proportional  to  the  number  of 
reversals. 


The  number  of  reversals  is  a meaningful  variable  as 
long  as  the  strain  amplitude  is  constant.  For  example, 
in  the  most  general  variable  strain  amplitude  situation, 
in  which  the  strain  amplitude  differs  from  reversal  to 
reversal,  the  number  of  reversals  at  each  strain  ampli- 
tude never  exceeds  one.  Thus,  number  of  reversals  los- 
es significance  as  an  influencing  variable.  Therefore,  to 
keep  the  formulation  simple,  the  rate  of  yield  strength 
change  (or  that  of  YSI)  will  be  assumed  to  be  propor- 
tional to  the  difference  between  the  saturation  and  cur- 
rent values  of  yield  strength  (or  that  of  YSI).  The  satu- 
ration YSI  will  be  assumed  to  be  a unique  function  of 
strain  amplitude  and  independent  of  prior  cyclic  his- 
tory. Furthermore,  the  coefficient  of  hardening/soften- 
ing (C ns ) will  be  assumed  to  be  a material  constant. 
Although  these  assumptions  will  lead  to  a somewhat  in- 
accurate prediction  of  the  transient  rates,  the  resulting 
model  will  be  simple,  and  it  will  always  approach  the 
correct  saturation  state. 

Dislocation  theories  and  observed  stored  energy 
changes  in  cold  work26'29  suggest  that  the  initial  part 
of  a reversal  after  prior  plastic  deformation  and  strain 
hardening  should  be  associated  with  some  recovery  (or 
slight  softening)  the  latter  part  of  the  reversal  where 

-6 R L.  Segall  and  ] M.  Finney,  “The  K alion  Between 
Physical  Properties  and  the  Obseived  Dislocation  Distribution 
in  Fatigued  Metals."  Acta  Metallurgica,  Vol  It  (July  196.1), 
pp  685-690 

2 7 1 ,.  M Clarebrough.  M.  V . Hargreaves,  and  M.  H.  l.oretto. 
Changes  in  Internal  h.'nergy  Associated  With  Recovery  and  Re 
crystallization  of  Metals  (Interscience.  196.1),  pp  6.1-1  21 . 

2®( , K 1 lalford . Stored  hnergy  of  Cold  Work  Changes  In- 
duced by  Cyclic  Deformation,  Ph  D.  Thesis  (University  of  IT 
linois,  1966). 

29A.  S.  Iyer  and  P.  Gordon,  "Note  on  the  Changes  in 
Stored  Energy  Produced  by  Reversals  of  Deformation,”  /ram 
actions  of  Metallurgical  Society  of  American  Institute  of  Min- 
ing. Metallurgical,  and  Petroleum  engineers.  Vol  52  (1952). 
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reversed  plastic  deformation  occurs  should  be  associa- 
ted with  work  hardening.  This  means  that  the  yield 
strength  (or  YSI1  would  initially  decrease  and  then  in- 
crease in  a reversal  as  reversed  plastic  deformation 
takes  place.  Adopting  this  view.  Figure  19  shows  a con- 
ceptual cyclic  variation  of  the  yield  strength  under  con- 
stant amplitude  cyclic  hardening,  cyclic  softening,  and 
cyclic  saturation  conditions.  To  use  a simple  quantita- 
tive characterization,  the  variation  of  YS1  within  each 
reversal  is  assumed  to  be  as  shown  in  Figure  20.  During 
the  initial  "linear”  elastic  portion,  the  YS1  is  assumed 
to  remain  unchanged.  At  the  end  of  the  “linear”  elas- 
tic part,  the  YS1  is  assumed  to  suddenly  reduce  by  a 
small  fraction  of  the  reversal.  The  fraction  Crc  by 
which  the  YS1  decreases  at  the  end  of  the  elastic 
part  is  designated  as  coefficient  of  cyclic  relaxation/ 
creep  and  is  assumed  to  be  a material  constant.  During 
the  latter  part  of  the  reversal,  the  YS1  is  assumed  to  in- 
crease linearly  with  stress.  The  value  of  YS1  at  the  end 
of  the  reversal  is  determined  by  the  model  for  revet  sal 
to  reversal  variation  in  YSI  discussed  in  the  previous 
paragraph.  The  resulting  expressions  for  the  YSI  at  any 
given  point  in  the  reversal  are  listed  in  Figure  20.  These 
represent  the  initial  YSIs  for  the  subsequent  reversal  if 
the  loading  direction  is  changed  at  the  current  position 
of  the  present  reversal.  The  stress-strain  path  in  each  re- 
versal is  completely  determined  by  the  initial  YSI  of 
that  reversal.  Thus,  a characterization  of  the  transient 
phenomena  requires  two  material  constants  (C’ns  and 
Crc),  the  determination  of  which  will  be  discussed  in 
Chapter  4. 

It  is  therefore  possible  to  quantitatively  characterize 
the  various  cycle-dependent  transient  phenomena  and 
the  observed  deviation  of  real  materials  from  Masing 
behavior  in  terms  of  a single  history-dependent  param- 
eter the  YSI.  This  approach  has  a reasonable  physical 
basis,  since  Ihe  YSI  represents  the  history -dependent 
variation  in  the  internal  friction  stress  of  the  material 
caused  by  microslruetural  changes  due  to  cyclic  plastic 
straining. 


A DETERMINATION  OF  MATERIAL 
^ PROPERTIES 

Depending  on  the  number  of  features  to  be  incor- 
porated, development  of  a cyclic  stress-strain  reponse 
model  as  discussed  in  the  previous  two  chapters  re- 
quires some  or  all  of  the  following  properties; 


1 . Monotonic  stress-strain  curve 

2.  Skeleton  stress-strain  curve 

3.  First  reversal  YR1  (or  YSI)  expressed  as  a func- 
tion of  monotonic  strain 

4.  Saturation  YRI  (or  YSI)  expressed  as  a function 
of  strain  amplitude 

5.  Coefficient  of  cyclic  hardening/softening 

6.  Coefficient  of  cyclic  relaxation/creep. 

The  first  property  the  monotonic  stress-strain  curve 
is  obtained  from  an  ordinary  tensile  test.  Since  nominal 
strains  exceeding  approximately  0.02  are  not  of  inter- 
est in  structural  analysis,  the  stress-strain  curve  ex- 
pressed in  engineering  units  (based  on  original  specimen 
dimensions)  is  adequate.  The  rest  of  the  properties  can 
be  termed  “cyclic”  properties  which  require  special  cy- 
clic tests. 

Cyclic  Tests 

Two  types  of  cyclic  tests  are  of  interest  in  determin- 
ing cyclic  properties: 

1 . Fully  reversed  constant  amplitude  tests 

2.  A multiblock  decremental  step  test. 

In  both  of  these  tests,  a uniaxial  specimen  is  subjected 
to  axial  strain  cycling  between  fully  reversed  specified 
limits,  and  the  asssociated  stress-strain  response  is  re- 
corded on  an  X-Y  recorder.  The  specimen  design,  its 
preparation,  and  the  procedures  for  testing  are  de- 
scribed in  the  forthcoming  ASTM  standards30  and  a 
currently  available  ASTM  publication.31  Therefore, 
such  details  are  omitted  here. 

In  the  case  of  fully  reversed  constant  strain  ampli- 
tude tests,  several  identical  specimens  are  strain-cycled, 
each  at  a different  strain  amplitude.  Although  it  is  suf- 
ficient to  cycle  the  specimens  until  saturation  is  reached, 
continued  cycling  to  fatique  fracture  yields  additional 
valuable  information  for  establishing  fatigue  properties 
Since  a real  test  may  never  reach  an  ideal  state  of  sat- 

3nRccommendcd  Practice  of  Constant  I mplitude  l ow  t i 
tie  Fatigue  Testing,  l-'orthcoming  ASTM  Standards.  Pari  In 
(ASTMl 

■'*  fl STM  Manual  on  Tow  Cycle  Fatigue  and  Testing,  ASTM 
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Soyi  + | = fiayj  + C|i  (6aysj  5ayi) 

Lei  Ao  = | a 0j| 

<S(jy  = 5a yj  forAa<  Aay 

5oy  = ( 1 Cr  ) 6ayj  for  Aa  = Aay 

5ay  - ( I cR)  5ayi  + [(Cr  cji)  5ayi  + C|i5aysi|  • \~  | for  Aa  > Aay 

where  Aay  = elastic  stress  range 

5oysj  = saturation  yield  strength  of  range  i I to  i 
c(|  = coefficient  of  cyclic  hardening  or  softening 

cr  = coefficient  of  cyclic  relaxation  or  creep 

5uy  = yield  strength  increment 


Figure  20.  A simple  model  for  the  cyclic  variation  of  yield  strength  (or  YSi ) 


oration,  such  a condition  may  he  defined  as  corres- 
ponding to  minimal  changes  in  the  stress-strain  hystere- 
sis loops.  In  some  cases,  the  stage  corresponding  to  one- 
hall  ol  the  fatigue  life  has  been  adopted  as  the  state  of 
saturation.-12.  11  Figure  21  shows  a typical  record  of 
stress-strain  hysteresis  response  obtained  for  A-36  steel 
during  the  first  four  cycles.  Since  changes  in  the  hyster- 
esis paths  are  significant  during  the  first  few  cycles,  a 

32R  W Lmdgrjf*.  Cvelie  Deformation  and  h'atigue  of  Har 
defied  Steels,  I'&AM  Report  .720  (University  oi  Illinois.  plfiKi 

33S  Keshas. in.  Some  Studies  on  the  Deformation  and 
I ras  lure  oj  V< mnaiitt'd  Steel  Under  Ceclic  Conditions  I'll  l> 
Dissertation  (University  o(  Waterloo,  IHveinlrer  l%hi 


continuous  record  is  obtained  during  these  cycles.  As 
cycling  continues,  it  is  adequate  to  record  hysteresis 
loops  at  logarithmic  increments  of  cycles  or  reversals. 

Figure  22  illustrates  a typical  strain  history  imposed 
m a multiblock  decremental  step  test.  A single  specimen 
is  strain-cycled  at  various  strain  amplitudes,  starting 
from  the  largest  magnitude.  At  each  strain  amplitude, 
the  cycling  is  continued  until  saturation  is  achieved  at 
that  level  and  the  saturated  hysteresis  loop  is  recorded 
on  an  \ Y recorder. The  next  lower  strain  amplitude  level 
is  attained  by  continuous  decrements,  as  illustrated  in 
figure  22  Thus,  it  is  possible  to  obtain  saturated  hys- 
teresis loops  corresponding  to  several  strain  amplitudes 
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Figure  21.  Typical  recor 
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from  a single  lest;  however,  this  tesl  does  not  provide 
any  information  on  the  transient  properties.  In  cases 
where  only  saturation  properties  are  needed,  this  test 
is  sufficient  to  derive  all  the  cyclic  material  properties. 

Determination  of  Cyclic  Properties 

Skeleton  Stress-Strain  Curve 

As  explained  in  Chapter  3,  the  skeleton  stress-strain 
curve  is  one-half  of  the  invariant  portion  of  a hysteresis 
branch.  Since  all  the  variation  caused  by  transient  and 
non-Masing  behavior  is  associated  with,  the  elastic  por- 
tion. the  skeleton  stress-strain  curve  is  comprised  of  the 
least  elastic  portion,  f'o  determine  (Ins.  it  w ill  be  neces- 
sary to  compare  all  the  Inst  reversal  hysteresis  branches 
m the  case  of  cyclic  hardening  materials,  such  as  2024- 
14  aluminum,  and  all  the  saturation  hysteresis  branches 
in  the  case  of  softening  materials,  such  as  cold-formed 
steels  In  the  case  of  materials  which  exhibit  either 
softening  or  hardening,  it  may  be  necessary  to  compare 
both  the  first  reversal  and  saturation  hysteresis  branches. 
However,  based  on  the  observation  of  three  normal 
sieels  SAL  1018  normalized  steel,  SAL  1045  normal- 
ized steel,  and  A-36  steel,  all  of  which  exhibit  either 
softening  or  hardening  it  was  sufficient  to  compare  the 
saturation  hysteresis  branches  of  these  materials  to  de- 
termine the  skeleton  stress-strain  curve. 

Figures  14  and  15  illustrate  a convenient  procedure 
for  determining  the  skeleton  of  stress-strain  curve  from 
a comparison  of  the  saturation  hysteresis  branches  of 
A-36  steel.  The  procedure  would  be  similar  in  the  case 
of  first  reversal  branches.  The  various  sized  hysteresis 
loops  are  superimposed  on  their  lower  (or  upper)  tips 
as  shown  in  Figure  14,  and  then  translated  along  the 
elastic  slope  until  all  the  upper  (or  lower)  branches  are 
matched  as  well  as  possible,  as  shown  in  Figure  15. 
Curve  OAB  defines  doubled  skeleton  stress-strain  curve. 

It  is  of  interest  to  note  the  difference  between  the 
skeleton  stress-strain  curve  and  the  cyclic  stress-strain 
curve  defined  in  low-cycle  fatigue  literature. 34  The  cy- 
clic stress-strain  curve  is  defined  as  the  locus  of  tips  of 
saturated  hysteresis  loops,  and  therefore,  the  curve 
joining  the  locus  ol  branch  tips  defines  the  doubled  cy- 
clic stress-siiam  curve.  Figure  23  illustrates  the  differ- 
ence between  the  skeleton  and  the  cyclic  stress-strain 
curves  for  the  three  types  of  materials.  The  difference 
in  stiess  levels  between  these  two  curves  at  a specified 
plastic  strain  level  is  the  saturation  YSI  corresponding 
to  that  plastic  strain. 

\\  I .milgral  . I Morrow,  amt  I Lndo,  "Determina- 
tion ol  Hie  ( yelk  Stress-Strain  Curve,"  Journal  oj  Materials, 
ISIMJMI.SA.  Vol4.No  | i March  1 969).  pp  1 76-1 88. 


When  approximating  the  skeleton  stress-stiam  curve 
by  piecewise  linear  segments,  it  is  desirable  to  choose 
arbitrary  lengths  for  segments  in  order  to  accurately 
represent  the  curve  with  a minimum  number  of  seg- 
ments; this  is  significant  for  reducing  computational 
cost.  Chapter  5 illustrates  two  examples  of  piecewise 
linear  approximation. 

First  Reversal  anil  Saturation 
YRIs  or  YSIs 

Since  the  skeleton  stress-strain  curve  is  associated 
with  the  least  elastic  portion,  by  definition,  the  YSI 
corresponding  to  the  skeleton  curve  is  zero:  hence,  YKI 
(or  YSI)  corresponding  to  any  other  hysteresis  branch 
is  always  positive.  I he  YR1  corresponding  to  any  given 
hysteresis  branch  can  be  determined  by  measuring  the 
difference  in  the  elastic  portions  between  the  hysteresis 
branch  and  doubled  skeleton  stress-strain  curve  in  stress 
units.  As  their  names  indicate,  the  first  reveisal  and  sat- 
uration YRIs  correspond  to  those  branches.  Since  each 
strain  amplitude  is  associated  with  a first  reversal  hys- 
teresis and  a saturation  hysteresis  loop,  the  first  rever- 
sal and  saturation  YRIs  can  be  expressed  as  functions 
of  strain  amplitude.  Figures  10.  II.  and  12  were  ob- 
tained in  this  manner. 

Coefficient  of  Cyclic  llardcnin.it/Softcninf ; 

As  discussed  in  Chapter  3.  the  coefficient  of  cy  clic 
hardening/softening  is  generally  not  only  a function  of 
the  material  but  is  also  dependent  on  strain  amplitude. 
However,  to  keep  the  calculation  simple,  this  coeffi- 
cient is  being  treated  as  a material  constant.  The  follow- 
ing procedure  may  be  used  to  select  a reasonable  value 
for  this  coefficient,  based  on  the  assumption  that  the 
rate  of  change  in  yield  strength  per  reversal  is  propor- 
tional to  the  difference  between  the  saturation  and  cur- 
rent values  of  yield  strength  presented  in  Chapter  3,  the 
following  expression  for  the  coefficient  of  cyclic  har- 
dening/softening, C|is.  can  be  derived. 


where  f = degree  of  saturation  expressed  as  a frac- 
tion (f  = 1 represents  complete  satura- 
tion) 

R = number  of  reversals  to  reach  this  degree 
of  saturation  under  constant  strain  cy- 
cling. 

For  example,  to  achieve  a 40  percent  saturation  in 
10  reversals,  a Cps  - 0.25  is  required.  It  is  possible  to 


pick  a reasonable  sel  of  I and  R values  from  the  ob- 
served transient  behavior  in  constant  strain  cycling 
tests.  In  the  absence  of  such  data,  a reasonable  value  ot 
(.'ns  based  on  the  observed  behavior  of  three  materials 
(see  Chapter  3)  appears  to  be  0.1  to  0.2. 

Coefficient  of  Cyclic  Relaxation /Creep 

As  explained  in  Chapter  3.  there  are  no  quantitative 
measurements  of  the  degree  of  recovery  on  unloading 
after  prior  plastic  deformation.  Thus,  knowledge  about 
the  coefficient  of  cyclic  relaxation/creep  related  to  this 
phenomenon  is  incomplete.  A value  of  0.1  to  0.5  ap- 
pears to  be  reasonable,  since  the  use  of  such  a value  in 
a model  results  in  a qualitatively  reasonable  simulation 
of  observed  cyclic  relaxation/creep  behavior  of  A-36 
steel  and  2024-T4  aluminum  (see  Chapter  5). 

Therefore,  the  various  material  parameters  (or  prop- 
erties) required  by  the  inelastic  response  model  can  be 
determined  from  a set  of  simple  cyclic  uniaxial  tests  of 
constant  strain  amplitude. 


C MODELING  CYCLIC  STRESS-STRAIN 
° RESPONSE 

I bis  chapter  discusses  three  aspects  of  the  model  for 
cyclic  stress-strain  response: 

1 . Development  of  a model  to  simulate  both  memory 
and  plastic  hysteresis  phenomena 

2.  Degrees  of  approximation/simplification  in  the 
model  and  appropriate  areas  of  application 

• 

3.  Typical  simulations  of  two  different  materials’ 
behavior. 

Combined  Model  for  Memory  and 
Plastic  Hysteresis  Phenomena 

The  memory  model  presented  in  Chapter  2 can  be 
simply  modified  to  incorporate  the  transient  and  non- 
Masittg  characteristics  of  real  materials  by  varying  the 
length  of  the  first  linear  segment  of  the  stress-strain 
curve  as  denoted  by  the  variation  in  YSI.  The  model 
foi  the  history -dependent  variation  of  YSI  was  dis- 
vti'ced  m I liaptei  3 I he  basic  length  and  slope  of  the 
liisi  segment  and  the  lengths  and  slopes  of  othei  seg- 
ments aie  defined  bv  the  skeleton  stress-strain  curve 
and  are  therefore  history-independent  material  con- 
stants 


cr  A E 


Figure  24.  Effect  of  relaxation  on  memory  phenomenon. 


In  addition  to  the  above  modification,  it  is  neces- 
sary to  modify  the  memory  rules  slightly  lor  the  availa- 
bilities of  certain  segments  during  relaxation/creep  to 
properly  simulate  observed  material  behavior.  Figure 
24  illustrates  this  problem  for  the  case  of  relaxation. 
The  difference  between  the  branches  AB  and  BC  is  due 
to  YRI  of  BC  being  less  than  YRI  of  AB.  II  loading  is 
continued  beyond  C,  the  memory  model  (or  rules)  will 
describe  a path  CD  parallel  to  AF,  which  is  a continua- 
tion of  OA;  however,  real  materials  follow  a path  simi- 
lar to  CD.35  Although  this  appears  to  be  a minor  dis- 
crepancy in  the  example  illustrated,  it  will  cause  signifi- 
cant errors  in  Ihe  simulation  of  cyclic  creep. 

The  observed  material  behavior  can  be  simulated  by 
recovering  the  availabilities  of  those  segments  that  are 
either  partially  available  or  totally  unavailable  to  the 
extent  of  decrease  in  YRI.  The  following  additional 
memory  rules  (memory  rule  5)  is  added  to  those  pre- 
sented in  Chapter  2. 


35ll  K Jliansale.  InelaMie  Deformation  an J langur  lie 
s/Hinse  of  Sprt  truni  I oaJeJ  Strain  ( ontrollea  t v talanJ  h'lexural 
Members,  Pill)  Diva-nation  (ttniversily  n!  W.ili-tloo,  I'll-I  I 


Cable  2 

Illustration  of  the  Application  of  Memory  Rule  5 to  the  Simulation  Example  for  Figure  27(a) 
for  the  Last  Reversal  Between  Strain  Limits  +0.01  and  +0.02  (A-36  Steel) 


Segment  No. 

t 

2 

3 

4 

5 

6 

7 

8 

Stress  increment  ot' 
segment  in  ksi  (MPa) 

20.00 
(137  90) 

10.00 

(68.95) 

5.00 

(34.47) 

5.00 

(34.47) 

5.00 

(34,47) 

2.50 

(17.24) 

1.00 

(6.89) 

1.50 

(10.34) 

Availability  ot  segment 
before  application  of 
memory  rule  5 

2.00 

2.00 

2.00 

2.00 

1.73 

0.00 

0.00 

0.95 

Maximum  stress  that  can 
be  recovered  in  ksi  (MPa) 

0.00 

(0.00) 

0.00 

(0.00) 

0.00 

(0.00) 

0.00 

(0.00) 

1.35 

(9.31) 

5.00 

(34.47) 

2.00 

(13.79) 

1.575 

(10.859) 

Actual  stress  that  is 
recovered  as  per  memory 
rule  5.  in  ksi  (MPa l 

0.00 

(0.00) 

0.00 

(0.00) 

0.00 

(0.00) 

0.00 

(0.00) 

0.49 

(3.38) 

0.00 

(0.00) 

0.00 

(0.00) 

0.00 

(0.00) 

Availability  of  segment 
after  application  of 
memory  rule  5 

2.00 

2.00 

2.00 

2.00 

1.83 

0.00 

0.00 

0.95 

Legend:  Given  A,.  It,.  and  21  l)j,  which  is  the  total  amount  of  relaxation  between  previous  and  current  reversals, 

the  quantities  Cj . l)j,  and  I j are  calculated  from  the  following  equations: 

1.  Cf  = (2.00  Bj)  Aj 

’ DjsC, 

8 

3 2.  D;  = 0.49  (Known) 

i = 1 1 

4 li=  13,+  Dj/Aj 


If  transient  features  are  simulated,  and  if  the  YRI  at 
the  beginning  of  the  current  reversal  is  less  than  that  at 
the  beginning  of  the  previous  reversal,  the  availability 
coefficients  of  the  segments  associated  with  the  current 
reversal  ate  increased  as  follows.  Starting  with  the  first, 
and  m consecutive  order,  the  availability  coefficients  of 
as  many  segments  as  required  are  increased  so  that  the 
total  sum  >)l  increases  in  stress  increments  is  equal  to 
the  difference  in  YRIs  between  the  beginning  of  the 
previous  and  current  reversals.  Hie  increase  in  uvaila- 
bility  coefficients  for  each  segment  is  limited  by  the 
condition  that  the  availability  coefficient  does  not  ex- 
ceed two  in  conformity  with  memory  rule  2.  The  avail- 
ability coefficients  of  segments  corresponding  to  oppo- 
site direction  are  appropriately  reduced.  The  applica- 
tion of  this  rule  is  illustrated  m Table  2 for  the  last  re- 
versal of  the  example  presented  in  Figure  27a,  which  is 


described  in  the  section  Typical  Models  and  Simulation 
for  Two  Materials. 

The  monotonic  stress-strain  curves  of  softening  and 
strain-aging  materials  are  different  in  shape  from  the 
rest  of  the  hysteresis  branches  described  by  the  skele- 
ton stress-strain  curve.  Therefore,  to  simulate  the  initial 
monotonic  path  accurately,  it  is  necessary  to  initially 
assign  the  various  segments  of  the  stress-strain  curve,  it 
is  then  necessary  to  change  these  values  to  those  given 
by  the  skeleton  stress-strain  curve  at  the  end  of  mono- 
tonic loading. 

Approximations  and  Appropriate 
Applications 

Whereas  a model  capable  ot  simulating  all  the  ob- 
served material  phenomena  is  generally  useful,  simplifi- 


Table  3 

I’ovdble  Degrees  of  Approxiinatioii/Simplification  in  Modeling 
Material  Behavior  and  Appropriate  Applications 


Model/ Degree  Material  Features  That 

of  Approximation  C an  Be  Simulateda 

Material  Properties  1 hat 
Are  Needed  and  the  Mini- 
mum Number  of  Material 
Characterization  Testsb 

Application s)  Which 
Require  This  Model 

1 All  four  features: 

A to  1) 

All  six  properties:  1 - 6 
Number  of  tests:  6 - 8 

Deformation  analyses  with 
relatively  few  cycles  of 
inelastic  loads 

2 l eatures:  B to  D 

Properties:  3 to  6 
Number  of  tests:  2 

Deformation  . na lyses  with 
several  cycles  of  loading 
Fatigue  analysis  where  mean 
stress  effect  is  important. 

3 l eatures:  C & 1) 

Properties:  5 & 6 
Number  of  tests.  1 - 2 

Deformation  analysis  where 
transient  behavior  is  unimpor- 
tant. Fatigue  analysis  involving 
non-Masing  materials  where 
mean  stress  effect  is  not 
important. 

4 Feature  D only 

Cyclic  stress-strain  curve 
Number  of  tests:  1 

Same  as  for  model  3.  except 
the  application  is  limited  to 
Masing  materials. 

“Materia/  Features:  A.  (mtial  cyclic  liardening/softening;  li  transient  behavior  due  to  irregular  cyclic  load- 
ing;/'. stable  non-Masing  behavior;  and  L).  stable  Masing  behavior. 

'’Material  Properties:  I Monotonic  stress-strain  curve;  2 first  reversal  of  yield  strength  increment; 3 coef- 
ficient of  cyclic  hardening/softening;  4.  coefficient  of  cyclic  relaxation/creep;  5.  skeleton  stress-strain 
curve;  and  6 saturation  yield  strength  increment  or  cyclic  stress-strain  curve. 


cations  or  approximations  may  still  provide  adequate 
models  for  certain  applications.  From  the  point  of  view 
of  such  simplifications,  the  observed  material  phenom- 
ena discussed  and  analyzed  in  the  previous  chapters  can 
be  divided  into  four  features: 

1 . Initial  cyclic  hardening  or  softening 

2.  Transient  behavior  due  to  irregular  cyclic  loading 
which  includes  combined  cyclic  hardeningor  soft- 
ening and  cyclic  relaxation  or  creep  (see  Chapter 

3) 

3.  Memory  phenomenon  with  non-Masing  charac- 
teristics 

4.  Memory  phenomenon  with  Masing  characteristics. 

Four  models  which  include  three  degrees  of  simplifi- 
cation result  from  the  above  categorization,  as  shown 
in  Table  3.  Model  I is  capable  of  simulating  all  the  four 


material  features  described  above  and  is  therefore  gen- 
erally useful  for  all  applications  listed  in  the  Iasi  column 
of  Table  3.  Models  2.  3,  and  4 correspond  to  increasing 
degrees  ot  approximation  or  simplification.  Model  2 
ignores  the  initial  cyclic  hardening  or  softening  feature, 
and  the  material  is  not  assumed  to  be  initially  at  the 
state  of  saturation.  Models  3 and  4 are  not  capable  of 
simulating  transient  behavior.  Although  mean  stresses 
are  important  in  some  situations,  transient  phenomena 
are  not  usually  considered  to  he  significant  in  cumula- 
tive fatigue  damage  analysis.  Therefore,  models  3 and  4 
would  be  adequate  for  use  in  fatigue  analysis  problems, 
where  mean  stress  effect  is  not  important. 

Table  3 also  lists  the  required  material  properties 
and  the  number  of  tests  to  determine  the  properties  of 
each  model.  Obviously,  model  1 requires  all  six  proper- 
ties (see  Chapter  4)  and  the  other  models  require  only 
some  of  these  properties.  It  is  of  interest  to  note  that 
the  number  of  material  characterization  tests  required 
to  formulate  model  2 are  about  the  same  as  those  for 


models  3 and  4,  Inti  model  2 is  much  superior  to  models 
3 or  4.  Also,  to  improve  from  model  2 to  model  1,  a 
considerable  increase  in  the  material  characterization 
tests  is  needed. 

Typical  Models  and  Simulation 
for  Two  Materials 

The  distinguishing  transient  behavior  features  of 
2024-14  aluminum  and  A-36  steel  will  be  illustrated  by 
appropriate  models  for  these  two  materials.  Figure  25 
illustrates  the  observed  cyclic  relaxation  behavior  of 
2024-T4  aluminum  and  a mild  steel  (which  is  similar  to 
A-36  steel)  under  strain  cycling  conditions  in  the  pres- 
ence of  a mean  stress.  The  2024-T4  aluminum  exhibits 
a simultaneous  cyclic  hardening,  whereas  tire  mild  steel 
(and  strain-aging  structural  steels  in  general)  exhibits  a 
simultaneous  cyclic  softening  in  addition  to  cyclic  re- 
laxation. These  phenomena  appear  as  cycle-dependent 
increases  or  decreases  in  stress  range,  respectively,  as 
shown  by  the  figure.  This  difference  can  be  attributed 
to  the  different  in  the  saturation  YSI  properties  be- 
tween the  two  materials,  this  will  be  demonstrated  by 
model  simulation. 

2024-T4  Aluminum 

Using  the  model  1 degree  of  sophistication.  Table  4 
lists  the  required  material  properties.  In  the  case  of 
202414  aluminum,  the  monotonic  and  skeleton  stress- 
strain  curves  are  identical.  The  first  reversal  YSI  and 
saturation  YSI  are  constants  and  independent  of  strain 
amplitude.  For  illustration,  the  skeleton  stress-strain  is 
approximated  by  nine  linear  segments;  the  coefficients 
of  cyclic  hardening/softening  and  coefficient  of  cyclic 
relaxation/creep  are  chosen  as  0.1 . 

Figure  26  shows  the  simulated  stress-strain  response 
of  2024-T4  aluminum  for  a specified  input  strain  his- 
tory. Figure  26a  illustrates  cyclic  hardening  during  the 

Ae 

first  eight  reversals  of  strain  amplitude,  — = 0.012. 

and  Figure  26b  illustrates  the  combined  cyclic  relaxa- 
tion and  hardening  in  the  presence  of  a tensile  mean 
stress. 

A-36  Sled 

Adopting  the  model  2 formulation  for  this  case,  the 
initial  cyclic  softening  or  hardening  leature  is  neglected; 
i.e..  the  material  is  assumed  to  be  in  the  saturated  con- 
dition initially.  The  initial  monotonic  stress-strain  path 
is,  therefore,  given  by  the  cyclic  stress-strain  curve.  Ta- 
ble 5 lists  the  material  properties  required  for  this  formu- 
lation Only  two  properties  the  skeleton  stress-strain 
curve  and  the  cyclic  stress-strain  curve  need  to  be  known. 


Table  4 

Cyclic  Material  Properties  of  2024-T4 
Aluminum  for  Model  1 


1 . Monotonic  stress-strain  curve  is  the  same  as  skeleton  stress- 

strain  curve. 

2.  Skeleton  stress-strain  curve  (piecewise  linear 

approxima- 

tion). 

Node 

Stress,  ksi  (MPa) 

Strain 

1 

0.0  (0,0) 

0.0 

2 

30.0  (206.9) 

0.003 

3 

35.0  (241.3) 

0.0037 

4 

411.0  (275.8) 

0.0050 

5 

42.5  (293.0) 

0.006 

6 

45.0  (310.3) 

0.0073 

7 

47.5  (327.5) 

0.0098 

8 

48.75  (336.13) 

o.oi  i 

9 

50.(1  (347.5) 

0.0131 

10 

51  2 (353.0) 

0.0153 

3.  1 irst  reversal  YSI 

= ().()  ksi  (().()  MPa). 

4 Saturation  YSI  = 

25.0  ksi  (172.4  MPa). 

5.  Coefficient  of  cyclic  hardening/softening  = 0.1. 

6.  Coefficient  of  cyclic  relaxation/creep  = 0.1 . 

Table  5 

Cyclic  Material  Properties  of  A-36  Steel  for  Model  2 

1 Skeleton  stress-strain  curve  (piecewise  linear  approximation) 

Node 

Stress,  ksi,  MPa 

Strain 

1 

0.0  (0.0) 

0.0 

2 

20.0  H 37.9) 

0.0007 

3 

30.0  (206.9) 

0.001  3 

4 

35.0  (241.3) 

0.0017 

5 

40.0  (275.8) 

0.0027 

6 

45.0  (310.3) 

0.0053 

7 

47.5  (327.5) 

0.0092 

8 

48.5  (334.4) 

0.015 

9 

50.0  (347.5) 

0.03 

2.  Saturation  YSI  (fioys)  in  ksi  expressed  as  a function  of  plas- 

tie  strain  amplitude  (Arp/2): 

A(  p Ai  p 

ftoys  = 1745  X /.lor  — -<  0.002. 

6oys  = 620  X A'1P  ♦ 2.25.  for  0.002  • ',P  • 0.0045. 

)><iy5  = 291.25  X A',P  + 3.73  tor  0.0045  • ^',P  < .011 
Soys  = 5.55  X ^ + 7.5  lor  -/C  .•  0.01 1 
3.  Coefficient  of  cyclic  hardening/sol  toning  0.1 
4 Coefficient  of  cyclic  relaxation/crcep  0.3 


(a)  Combined  cyclic  relaxation  and  softening  of  mild  (b)  Combined  cyclic  relaxation  and  hardening  of  2024 
steel  under  constant  strain  cycling.  T4  aluminum. 

Figure  25.  Cyclic  relaxation  behavior  of  2024-T4  aluminum  and  mild  steel. 


(a)  Cyclic  hardening.  (b)  Combined  cyclic  relaxation  and  hardening. 

Figure  26.  Simulation  ol  cyclic  hardening  and  combined  cyclic  relaxation  and  hardening  of  2024-T4  aluminum  un- 
der strain  cycling  history:  ( I ) 0 to  0.01 2,  one  to  eight  reversals  at  +0.01  2 and  (2)  nine  to  ten  reversals  at 
'0.01  2.  one  to  nine  reversals  at  +0  012  and  0.002  and  one  reversal  at  0.002  and  +0.01  5 


since  the  third  property  (saturation  YSI)  can  be  deter- 
mined front  the  oilier  two.  For  illustration,  both  the 
skeleton  and  the  cyclic  stress-strain  curves  are  approx- 
imated by  eight  linear  segments;  the  coefficient  ofhard- 
ening/softening  and  coefficient  ol  relaxation/creep  are 
chosen  as  0.1  and  0.3.  respectively. 

Figure  27  illustrates  the  simulated  cyclic  relaxation 
and  cyclic  creep  behavior  under  strain  and  stress  cycling 
histories,  respectively.  In  both  cases  the  model  exhibits 
a cyclic  softening  simultaneous  with  relaxation  or  creep, 
which  is  consistent  with  the  observed  behavior  of  mild 
steels. 

Figure  2N  compares  the  simulated  and  experimental- 
ly observed  stress-strain  response  under  a variable  strain 
history;  the  agreement  is  reasonably  good. 


£ APPLICATIONS  AND  SCOPE 
V FOR  FUTURE  WORK 

The  work  presented  in  this  report  is  directly  or  indi- 
rectly useful  in  three  applications  of  structural  mechan- 
ics and  failure  analysis. 

Basic  Material  Parameters 
for  Constitutive  Equations 

The  cyclic  material  properties  developed  in  this  work 
provide  the  necessary  basic  properties  for  formulating 
material  constitutive  equations  for  the  mulliaxial  stress 
states.  Such  constitutive  equations  are  directly  adopta- 
ble  in  finite  element  computer  codes  for  nonlinear 
structural  analysis.  As  explained  in  Chapter  1 . the  need 
lor  nonlinear  structural  analysis  techniques  in  seismic 
design  is  well  recognized.  A plasticity  formulation  for 
developing  appropriate  constitutive  equations  for  cyclic 
loading  conditions  which  uses  the  material  parameters 
developed  in  this  report  has  been  proposed. 3(1 

Direct  Hysteresis  Model 

for  Structural  Compouent/System 

A material's  inelastic  stress-strain  response  has  a 
major  influence  on  a structural  component's  load- 


36S  h Slurtna  and  II  K Jhansale.  I Plasticity  formula- 
non  for  Cyclic  Inelastic  Structural  \ Italy  sis.  Interim  Report 
M-202/APA037U47  (l  S Army  Construction  Engineering 
Research  laboratory.  1977). 


deformation  characteristic,  in  addition  to  other  factors 
such  as  stress  or  strain  gradients,  geometry,  or  joint  slip. 
The  hysteretie  load-deformation  behavior  often  has 
features  qualitatively  simiiar  to  the  material  stress-strain 
characteristics  discussed  in  this  report.  Therefore,  in 
such  cases,  cyclic  inelastic  load-deformation  models  can 
be  directly  formulated  along  the  lines  ol  development 
for  the  material  model  outlined  in  this  report.  The 
terms  "stress."  “strain,"  and  “cyclic  material  param- 
eters" in  the  material  model  are  analogous  to  the  terms 
"load,"  "deformation,"  and  “cyclic  components/system 
parameters”  in  the  structural  response  model.  The 
component  parameters  (or  properties)  required  by  the 
structural  model  can  be  determined  from  simple  cyclic- 
tests  of  components  similar  to  those  outlined  in  this 
report. 

A structural  reponse  model  would  be  directly  usetul 
in  the  area  of  seismic  analysis;  a few  applications  are: 

( 1 ) development  of  accurate  inelastic  response  spectra, 

(2)  energy  dissipation  due  to  inelastic  hysteresis,  and 

(3)  evaluation  of  realistic  ductility  factors. 

Fatigue  Analysis  Under  Complex  Histories 

Extensive  low-cycle  fatigue  research  during  the  past 
two  decades  has  shown  that  the  fatigue  process  consists 
of  an  initial  crack  initiation  stage  and  a final  crack 
propagation  stage.  The  stress  and  plastic  strain  fields  in 
a structure's  critical  locations  govern  the  rate  of 'damage 
during  these  two  stages  of  fatigue.  Therefore,  the  ma- 
terial model  discussed  in  this  report  is  directly  useful 
for  evaluating  the  accumulation  of  fatigue  damage  at 
specified  critical  locations  of  structures,  because  it 
provides  complete  information  on  the  stress-strain  re- 
sponse of  the  material  at  these  locations.  However,  the 
inelastic  response  model  plays  a more  important  role  in 
accounting  for  the  "sequence  effect”  in  fatigue,  which 
is  essentially  a manifestation  of  the  prior  history  depen- 
dence of  the  rate  of  damage  accumulation.37  Simple 
linear  damage  summation  procedures  which  ignore  this 
effect  have  resulted  in  extremely  nonconservative  or 
conservative  estimates  of  fatigue  lives.  Thus,  history- 
dependent  inelastic  response  models  would  be  extreme- 
ly valuable  for  achieving  more  accurate  fatigue  life 
prediction  procedures. 


37||.  it.  Jhansale,  “Evaluation  ol  Deformation  Phenomena 
of  Metals  for  Fatigue  Analysis."  Journal  of  Tcstiny  and  I va/tc 
alton.  ASTM-JTCIA,  Vol  3.  No.  5 (ASIM.  September  19751. 
pp  348-354. 
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(a)  Combined  cyclic  relaxation  and  softening  under  (b)  Combined  cyclic  creep  and  softening  under  stress 
strain  cycling  between +0.01  to +0.016.  cycling  between  -32.0  to  +54.00  ksi  (-220.6  to 

+372.3  MPa). 

Figure  27.  Simulation  of  transient  behavior  of  A-36  steel. 
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Figure  2N  Comparison  ol  simulated  and  experimentally  observed  hysteresis  behavior  of  A-36  steel  under  a variahl 
strain  history. 
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Conclusions 

The  observed  history-dependent  inelastic  behavior 
ot  structural  metals  having  a wide  range  of  strengths 
I and  cyclic  characteristics  can  be  accurately  modeled  by 

slightly  modifying  a spring-slider  rheological  model 
(f  igure  3)  or  its  equivalent  memory  model. 

The  various  cycle-dependent  transient  phenomena, 
including  hardening,  softening,  relaxation,  and  creep, 
and  the  observed  deviation  from  Musing  behavior  can 
be  quantitatively  characterized  in  terms  of  a single  his- 
tory-dependent yield  strength  parameter.  Although  this 
development  is  essentially  based  on  phenomenological 
observations  of  a large  number  of  structural  metals,  it 
appears  to  be  physically  consistent  with  mierostruetural 
processes. 

All  the  material  parameters  required  by  the  model 
can  be  determined  from  simple  cyclic  tests  of  constant 
• strain  amplitude. 

j i 

Impact  of  Study 

This  research  program  has  achieved  an  important 
step  in  formulating  a nonlinear  inelastic  structural  anal- 
ysis procedure  for  such  complex  loadings  as  those  ex- 
perienced during  earthquakes.  Incorporating  the  devel- 
opments of  this  research  should  make  it  possible  to  ac- 
curately predict  the  nonlinear  dynamic  response  of 
seismic  structures,  in  addition,  it  should  now  be  possi- 
ble to  evaluate  damping  coefficients  and/or  ductility 
factors  used  in  seismic  design  (see,  for  example,  TM  5- 
809-10.  Seismic  Design  for  Buildings ) on  a more  ra- 
tional basis  than  using  past  experience.  Thus,  the  study 
has  contributed  to  improvement  of  seismic  analysis  and 
design  methods  both  in  terms  of  reliability  and  cost  ef- 
fectiveness. Extension  of  this  research  into  fatigue, 
fracture,  and  shakedown  problems  would  provide  accu- 
rate quantitative  procedures  for  evaluating  time-lustoi  v- 
dependent  cumulative  damage  criteria  mulct  seismic 
loadings. 
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